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Abstract The Newton Modified Barrier Method (NMBM) is 
applied to structural optimization problems with large a number 
of design variables and constraints. This nonlinear mathematical 
programming algorithm was based on the Modified Barrier Func- 
tion (MBF) theory and the Newton method for unconstrained 
optimization. The distinctive feature of the NMBM method is 
the rate of convergence that is due to the fact that the design 
remains in the Newton area after each Lagrange multiplier up- 
date. This convergence characteristic is illustrated by application 
to structural problems with a varying number of design variables 
and constraints. The results are compared with those obtained by 
optimality criteria (OC) methods and by the ASTROS program. 


1 Introduction 

The problems facing designers today tend to be complex and 
often novel with limited past experience to guide the designer 
to satisfy conflicting requirements. The computer automated 
analysis and simulation capabilities have reached a very high 
level of sophistication and power to provide detailed infor- 
mation about the response behaviour of the design. These 
capabilities, however, often provide little information as to 
how to change a complex design to improve its behaviour or 
its cost measure, or both. New automated optimization ap- 
proaches are needed to provide guidance for improving the 
preliminary and for fine tuning the final designs. 

The use of finite element methods in structural analysis 
leads to a large number of size variables. It is often conve- 
nient to treat each element size as an independent variable 
instead of developing a variable linking procedure to reduce 
the number of variables, or to promote smoothness of member 
size variations. Such reductions of problem size are, however, 
always advantageous from a computational point of view. 

The derivations of OC methods (see Berke and Khot 1987 ; 
Khot 1981) are associated with structural energy theorems 
that are involved in stiffness related optimization problems. 
These energy theorems assure separability and the associated 


simplifications. The number of needed iterations depends not 
on the number of design variables but on low sensitivity of 
internal forces to changes in member sizes. These OC meth- 
ods have been shown to be extremely efficient but only for 
special cases of structural constraints. With the increasing 
importance of multidisciplinary design requirements special- 
ized algorithms are becoming less attractive for the emerging 
complexities. 

With the increasing speed of computers, a major part of 
the importance of computational efficiency needed for opti- 
mization iterations is diminishing. In such a computational 
environment, developments turn towards robust optimization 
methods based on sound mathematical foundations. 

The NMBM method applied here initially only to struc- 
tural problems promises to satisfy these requirements. Its 
performance is compared here with that of OC and the AS- 
TROS programs for problems that satisfy the conditions of 
the applicability of OC methods. Separability and the lin- 
ear objective function are not necessary conditions for the 
proposed method, but are believed to improve convergence 
behaviour. Work is in progress to compare this approach to 
other algorithms in the multidisciplinary problem setting us- 
ing COMETBOARD, a testbed for optimization algorithms 
(Gendy et ai 1994). 

The NMBM is a numerical realization of the general 
modified barrier method (MBM) for constrained optimiza- 
tion (Polyak 1992a). The barrier function (MBF) used in 
this method possesses the best properties of the classical La- 
grangian and classical barrier functions (CBF) (see Fiacco 
and McCormick 1968), and at the same time eliminates their 
main drawbacks. The MBF retains the best properties of the 
augmented Lagrangian (see Powell 1969; Rockafellar 1973; 
Bertsekas 1982) and also develops new important qualities. 
The properties of MBF near the primal and dual solution 
make it possible to avoid illconditioning of the CBF s Hes- 
sian. This makes the numerical process more stable and im- 
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proves the rate of convergence. The application of MBF-type 
methods for solving truss topology design problems (Ben-Tal 
et al. 1992) leads us to believe that it has good potential for 
also solving large-scale structural optimization problems. 

The next sections give the formulation of the structural 
optimization problem using the NMBM approach, followed 
by presentations of MBF, MBM, NMBM in turn. In the last 
section, three truss structures using NMBM are presented. 
The convergence characteristics are compared with solutions 
obtained by ASTROS (Neill and Herendeen 1993) and by an 
OC method. 

2 Structural optimization problem 

The structural optimization problem, which we are going to 
consider, can be stated as follows: 
n 

W* = min W = £ Pi t iX{ , (1) 

i=l 

subject to constraints 
n 

9j{x)=Cj-Cj = Y^QijX~ l -Cj< 0, j = 1 m, (2) 

i=l 

and 

x,* > 0, i = 1,... ,m. (3) 

The are the design variables, p t - is the specific weight, and 
is a function of the geometry of the elements. Parame- 
ters depend on the nature of the constraints and are a 
function of the sensitivity, and cj and Cj are the actual and 
limited values of the ;-th constraints. 

Before we start our analysis, let us rewrite the problem 
(l)-(3) in a more convenient manner for further consideration. 
We introduce the reciprocal variables y^ = x“*,t = 
and the vector y = (y/,...,yn)* We also set 
r, = so that the objective function (1) can be rewritten 
as follows: 

w = T, r iyf 1 ’ 

i=l 

and for the constraints (2) we have 
n 

9j{y) = Cj~Yl Qj' y i - °> j - 1 m 

1=1 

The objective function 

W = f2r i y~ 1 
1=1 

is a barrier type function; therefore, if the initial approxima- 
tion y° = (yj , . . . , y£) is inside the feasible set, i.e. ail y? > 0, 
then the minimization procedure will keep the variables pos- 
itive. 

Let Q =|| ||, j = 1, . . . , m, i = l,...,n, be an 

m x n matrix with elements Qj t \r = (17, . . . , r n ) and y^ 1 = 

(yj” 1 j ■ * ■ 1 yn l ) ^e n-dimensional vectors, A = (A] , . . . , A m ) 
an m-dimensional vector of Lagrange multipliers, and g(y) = 

( z \ ~ E"=i QiiVi . ■ ■ ■ - £"=1 CmiVi) = c - Qy an Tri- 

dimensional vector function. Then the problem (l)-(3) can 
be written as follows: 

n 

W* = f(y*) = min/(y) = min(r,y _1 ) = ^r-,2/' 1 , (4) 

«=1 


subject to the constraints 

g(y) = c - Qy > 0 . (5) 

For structural design problems, the feasible set that is defined 
by (5), namely, 

n = {y : 9,{ y) > 0, j = , 

has a nonempty interior; i.e. there exists y° > 0 such that 

?/(y°)>0. j= 1 m. (6) 

Therefore, the Karush-Kuhn-Tucker (KKT) conditions hold 
true. 

Let R = [diag r,]^ be a diagonal matrix with entries 
r f *, x = and let L(y, A) be the Lagrangian for the 

problem (4)-(5) 

m 

£(y.A) = /( y) - x )9](y ) = (r.y -1 ) - (A,c-Qy) , 

J=1 

then, for any KKT pair (y,A), we have, 

V y L(y, A) = -Ry- 2 +AQ = 0, (7) 

c- Qy > 0,(A,C - Qy) = 0, A > 0. (8) 

The modified barrier method, which we describe below, gen- 
erates primal and dual sequences that converge to a KKT 
pair (y*, A*). 

Due to the nature of the structural design problem, the 
feasible set J? is bounded. Therefore, the objective function 
/( y) is strongly convex on Q. Thus the solution y * of the 
problem (4)-(5) is unique, while the set of Lagrange multipli- 
ers A that satisfies (7) and (8) is bounded due to the Slater 
condition (6). 

In the next section, we will introduce the modified barrier 
function for the structural optimization problem (4)-(5). 

3 Modified barrier function 

Let k > 0 be a scaling parameter. The constraints (5) can be 
rewritten in the following form: 

k~ l ln[*yy(y) + 1] > 0, (9) 

We shall assume that In t = — oo for t < 0. 

The modified barrier function is a classical Lagrangian for 
the equivalent problem (4) and (9) 
n m 

F(y,\,k) = ^r.y" 1 - fc -1 ^ A j ln [*0y(y) + l l- ( 10 ) 

‘=1 J=1 

The MBF theory and the corresponding methods for nonlin- 
ear and linear optimization have been developed by Polyak 
(1992a, b). For the sake of completeness, we shall specify the 
main MBF properties for the structural optimization problem 
(4)-(5). 

Let Y = (diag be a diagonal matrix in n- 

dimensional space with entries y,* and A = (diag AyJJJLj be a 
diagonal matrix in m-dimensional space with entries A j. 

For any KKT pair (y* r A*) that satisfies (7) and (8), and 
any Jk > 0, the following properties are true: 

n 

PI • F(y*, A*, k) — /o(y*) = . 

1 = 1 

P2 *V y F(y*, A*, k) = V y L(y*,A*) = -R(y*)" 2 +A*Q = 0 
P3 • V 2 y F(y*,A*,t) = 2R(Y*) -3 + jfcQ T ^*Q , 


I 

i 

I 

* 


* 

* 

% 


« 
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where (y*)~ 2 = [(yp -2 , . . . .(yj)" 2 ] is a vector in R n and 
(Y*)~ 3 = [diag(y*)“ 3 ] [L.J is a diagonal matrix with entries 

(y*r 3 - 

The properties P1-P3 show that in contrast to CBF, the 
MBF as well as its first and second derivatives exist at the pri- 
mal solution and possess important properties in the neigh- 
bourhood of the primal-dual solution. 

The optimal value of the objective function coincides with 
the MBF value for any KKT pair (y*, A*) and y* is a unique 
minimizer of F( y, A*, k) for any k > 0. 

Taking into account that during the minimization proce- 
dure the vector y is positive, the MBF F(y, A.jfc) is strongly 
convex and smooth in an extended feasible set, = {y : 
9j( y) > 1, . . . ,m}. 

The MBF properties enable us to develop a method, 
which converges to the primal-dual solution for any positive 
parameter k > 0 (see Jensen and Polyak 1994). This means 
that in contrast to the interior point methods, that are based 
on classical barrier functions (see Gonzaga 1992), the MBF 
approach provides convergence without increasing the barrier 
parameter unboundedly. 

Therefore, the condition number of the MBF Hessian is 
stable from the beginning to the end of the process. It con- 
tributes to both the rate of convergence and numerical sta- 
bility. It also makes more efficient use of Newton's method 
for unconstrained minimization F(y, \,k) in y, which is the 
key element in the numerical realization of the MBF method 
that we are going to describe in the next section. 

4 Modified barrier methods 

The MBF method consists of unconstrained^ minimization 
F(y, A,*) in y and updating the vector of Lagrange mul- 
tipliers A. The barrier parameter k > p can be fixed or one 
can change it after every Lagrange multiplier update. We 
will now describe two versions of the MBF method, starting 
with the permanent parameter version. 

We start with positive vector y 3 , which belongs to the 
extended feasible set, i.e. 

9j( y°)>-* -1 , j = (11) 

It is worth mentioning that for any positive vector y°, one 
can find such k > 0 that inequalities (11) will be satisfied. 

Thus, finding the initial primal approximation does not 
require extra computational work. As far as the initial vec- 
tor of Lagrange multipliers is concerned, we can choose any 
positive vector A 0 = (Aj , . . . , A^). 

Let k > 0 be fixed, y £ 1?*, and A 0 = (1, . . . , 1) £ R n . 
Assume that the approximations y s and X s have been found 
already. Then the next approximations for y J+1 and A 5+1 
are found by the formulae 

y ,+1 :F(y s+t ,X s ,k) = mm{F(y,\ 3 ,k) |yeR n ), (12) 
and 

\ s+1 ■.y +l = \ a J [kg J (y’+') + l)-\ j = 1 , . . . , m . (13) 

From (11) and (12), we have 

Vy/V +1 , X s , k) = V y I(y J+1 , A' +1 ) = 0 , (14) 

i.e. the unconstrained minimum of the MBF F( y, A 3 ,Jfc) co- 
incides with the unconstrained minimum of the classical La- 
grangian L(y, A J+1 ) in y. 


It was proven (see Jensen and Polyak 1992) that if 17 
is bounded and (6) holds true, then the method (12)-(13) 
is executable, i.e. the unconstrained minimizer y s exists for 
any k > 0. In the case of the structural optimization problem 
(4)-(5), the primal minimizer y 3 is always unique. 

It was also proven that the method ( 12)-( 13) converges, 
i.e. lim 5 — oo/(y J ) = /( y*) whether or not the objective 
functions and the constraint functions are strongly convex. 
Moreover, the next estimation, 

l/(y s )-/(y*)l=(*a,r 1 /2 i (15) 

holds true for any fixed k > 0 and limj— oo <*3 = 0- For 
a structural design problem, the estimation (15) can be im- 
proved. The objective function f(y) for the structural design 
problem (4)-(5) is strongly convex; therefore, if A* is unique, 
then (due to Polyak 1992a) the MBF method ( 12)-( 13) con- 
verges with a linear rate. In other words, for any fixed and 
large enough k > 0 there exists 0 < q < 1 such that the 
following estimation holds true: 

max{|| y' +1 -y* ||,|| A ,+1 - A* ||} < 

(c/k) || A 1 - A* ||= q || A* - A* || , (16) 

and q — ► 0 if k — ► oo. Therefore by increasing the barrier 
parameter k step by step, one can improve the rate of con- 
vergence. 

Let us consider the second version of the MBF method 
y* +1 :F(y s+l ,X s ,k) = mi n {F(y,X s ,k s )\yeR n }, (17) 
A ,+ 1 : A* +1 = Aj[* aS y(y ,+ 1 ) + 1] _1 , j = l, ... ,m . (18) 
If A* is unique, then the following estimation is true: 
max{|| y 3 - y* ||,|| X 3 - X* ||} < 9l ■ . . . • q s , 
with q 3 — ► 0 if k 3 — ► oo. 

Before we turn to the numerical realization of the MBF 
method, we would like to emphasize that for every s > 1, 

V y L(y 4 5 ,A J ) = 0. (19) 

Also, 0,'(y J ) — *■ 0 for the active constraints, i.e. for j £ J = 

{j ■ s/(y*) = o}. 

For the passive constraints j : gj(y*) > <t > 0, therefore 
due to the formula for the Lagrange multipliers update, we 
have 

A* = A'-^y*) + H" 1 < X 3 (kff + I)" 1 , (20) 

i.e. A* < Xj(k<r + 1) — J — * 0,jcJ. 

We can see that, in contrast to the augmented Lagrangian 
method (see Powell 1969; Rockafeliar 1973; Bertsekas 1982), 
the MBF does not require any special care to keep the La- 
grange multipliers non-negative. It happens automatically 
because of (13) or (18) for Lagrange multiplier update. 

Due to (14) and the formula for Lagrange multiplier up- 
date, the sequence {(y J ,A J )}^ 0 converges to a KKT pair 
(y*,A*). To estimate the “distance” between the current 
approximation (y, A) and the solution (y*, A*) we introduce 
the following merit function: 

v(w,fc) - u(y, A, k) = 

m 

max{G(y), || V y L(y, A) ||, £ A,- | 9j {y) |} , 

;=1 

where G(y) = -{minify) | 1 < j << m}. 
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It is easy to see that for any feasible vector y. any positive 
Lagrange multiplier vector A = ( A i , . . . , A m ) and any k > 0 
the merit function is non-negative. 

It is also clear that 

v(y,A,*) = 0 (21) 

is equivalent to (7), so that 

i;(y,A,*) = 0 ^ ► y = y*, A = A* , 

for any k > 0. 

Also, there exists M > 0 such that 
v{s(w 3 f k) = v{y s ,\ 3 ,k) = 

V(y\ A 3 , k) - v(y\ \\ k) < M II A 5 - A* || . 

Therefore for any given 0 < 7 < 1 there exists k 1 > 0 such 
that 

v(w 3 ,k) = v(y 3 ,A 3 ,*) < y s , (22) 

for any k > k 7 . 

We shall use the merit function v(y , A , k) and the estima- 
tion (22) for the numerical realization of the MBF method. 

5 Newton modified barrier method 

The modified barrier methods (12)-(13) or (17)-(18) require 
solving the unconstrained optimization problem (12) or (17) 
at every step, which in turnj requires an infinite number of 
arithmetic operations. Therefore, to realize the MBF method 
numerically, we must find a way to replace the vector y 3 in 
(13) by such an approximation y 3 , which can be found in fi- 
nite number of arithmetic operations and retains the estima- 
tion (16). Such a way has been indicated by Polyak (1992a, 
see p. 205). 

In this paper, to find an approximation for y 3 we will use 
the Newton method up to the point when the norm in the 
Newton direction becomes smaller than a given small enough 
number e > 0, which is defined by the accuracy which we 
finally want to achieve. 

The Newton method is particularly efficient for structural 
optimization due to the following three facts. 

1. The MBF F(y, A, k) is smooth and strongly convex for y £ 

n k . 

2. It is easy to compute the MBF Hessian V yy F{y, A, k). 

3. The numerical linear algebra for solving normal systems of 
equations to find the Newton directions is well-developed. 

To describe the Newton MB method we introduce the 
matrix A(y, k) = {diag(ty J (y)+l]}''l 1 with entries kgj(y)+ 
1. Let us also remember that 

R = [diag rj]£L,,Y = [diag y,]" =1 ,A= [diag A,]”!, . 

For the gradient, we obtain 

VyF(y, \,k) = -Ry -2 - A^ _1 (y, ik)Q , (23) 

while the Hessian is given by 

V 2 y F(y, A, k) = 2Ry~ 3 + kQ T AA~ 2 (y , *)Q . (24) 

Now we will describe the NMBM. The NMBM has two basic 
parts. First, we find the approximation y 5 "^ for the primal 
minimizer y 3+1 . To find the approximation we will use the 
Newton method for minimizing the function F(y,A 3 ,Jb) in 
y when both the vector of Lagrange multipliers A 3 and the 
parameter k > 0 are fixed. 

Second, we update the Lagrange multipliers by using the 
approximation y 3+1 instead ofy 3+i in (13). 


After every Lagrange multiplier update, we check the 
merit function v(w,k), which is supposed to converge to zero 
at an a priori given rate 0 < 7 < 1. If the barrier parameter 
k > 0 is chosen properly, i.e. k > k 7 > 0, then every La- 
grange multiplier update will decrease the function v(w,k) 
by a given factor 0 < 7 < 1, i.e. 

v(w s+1 ,k)<v(w 3 ,k)y. (25) 

If the inequality (25) does not hold, it means that the barrier 
parameter is not large enough and we must increase it. 

The most costly operation in the NMBM is solving the 
system of equations 

V 2 y F(y,X 3 ,k)d=-V y F(y,X a ,k) l (26) 

to find the Newton direction d = d(y,A 3 ,fc). To guarantee 
convergence from any starting point y° £ J?, we must use 
the Newton method with a step size. 

To find a proper step size, we will use the relaxation in- 
equality, 


F(y + td, A, k) - F( y, A, k ) < (l/3)< [V y F(y, A, k), d] , (27) 


starting with t = 1. 

If (27) is fulfilled for t = 1, then we have a “pure” Newton 
step; otherwise, we check (27) for t := (i/2) and continue to 
do so while the inequality (27) will be fulfilled. 

Due to the properties of the MBF F( y, A, k), for any fixed . 

positive vector A > 0 and Jb > 0, it will come to a point where * 

the inequality (27) will be fulfilled for t = 1. At this point, we 
have reached the Newton area and, from this point on, we will 
use only “pure” Newton steps. Therefore, the minimization 
procedure will converge to the primal minimizer 

y(A, k) = argmin {^(y, A, k) \ x e R n } 

quadratically. The stopping criterion is given by the formula * 

II d HvJ y F(yA*) = H d HvJyF(-)— £ - ( 28 > i 


where c > 0 is the small enough number defined by the 
accuracy of the approximation to the solution y* we want 
to achieve. Now, we describe the Newton modified barrier 
method (see flow chart). 

Let s be the accuracy which we want to achieve for the fi- 
nal approximation y to the solution y*. We consider a mono- 
tonically increasing sequence : — +°°* 

We start with y° £ f?jfc,A 0 = (1 ,..., 1) £ R m ,jb = 
k( 0) = ko,l(Q) = 1, and 0 < 7 < 0.5. Now assume that 
y 3 , A 3 , fc(s),£(s) have been found already. To find the next 
approximation, one must perform the following steps. 

0. Start with y := y 3 . 

1. Set A := A 3 ,ib := k(s),£ := l(s). 

2. Find the Newton direction d = d(y,A,Ar) by solving the 
system 

V 2 yyF(y,\,k)d = — VyF(y,A,fc) (29) 




i 

i 

1 

1 



and set t = I. 

3. If y + td £ Q k , then go to 4; if y + td £ then t := (i/2) 
and go to 3. 

4. Check the inequalities (27). If they are not fulfilled, set 
t := (i/2) and check (27) again; if (27) is fulfilled and 
i = 1, go to 5; if (27) is fulfilled but 0 < t < 1, then set 
y := y + td and go to 2. 

5. If || d €, go to 6; otherwise, go to 2. 
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Fig. 1. Flow diagram for the NMBM method 


6. Set y := y,A = A— l(y,Jb)A,w = (y,A) and check the 
inequality 

< c . (30) 

If (30) is fulfilled, then w* = w and quit; else, go to 7. 

7. Check the inequality 

*(w)<7* +1 - (31) 

If (31) is satisfied, then set y 5+1 = y and A ,+ 1 = A, k(s + 
1) = jb(s),f(s + 1) = £(s) 4* l,s 4* 1 := s, and go to 1. If 
(31) is not satisfied, then set y 5+1 = argmin {/(y 1 ) | i = 
l,...,s+l},A i+1 = e = (l,...,l),fc(s+l) = fc 3 + 1 ,f(s + 
1) = £($), s 4- 1 := s, and go to 1. 

For any 0 < y < 1 there is a number s$ such that for s > 
sq, the barrier parameter will not change and the step-size 
t = 1; i.e. we will perform only the “pure” Newton method 
and the primal minimizer will stay in the Newton area after 
every Lagrange multiplier update. Such a point we call a 
hot-stari . 


From this point on, every Lagrange multiplier update 
requires at most 0(lnlne” 1 ) Newton steps, and every La- 
grange multiplier update shrinks the distance between the 
current approximation and the primal-dual solution by a fac- 
tor 0 < 7 < 1. This follows from the basic properties of MBF 
(see Polyak 1992a) and Newton’s method (see Smale 1986). 

The number sq depends on the condition of the con- 
strained optimization problem (see Polyak 1992a, b), which 
in turn depends on the “measure” of nondegeneracy of the 
problem (4)-(5). The value sq can be decreased by increas- 
ing the barrier parameter k > 0. Moreover, for any non- 
degenerate constrained optimization problem and any fixed 
0 < 7 < 0.5, there exists k y > 0 such that for any k > b 7 , 
one can achieve the hot-stari with sq = 1. 

However, if Jb > 0 is large enough it will take extra compu- 
tational work to reach the hot-start at the first stage, where 
the Lagrange multipliers are fixed and one must change the 
barrier parameter k > 0. To reach the hot-start , one can use 

shifted barrier function 
n 

F(x, e, k) = Yi W-' + 51 + J 1 • 

»=1 

In the case of the structural optimization problem, the 
shifted barrier functions possess the so-called self-concordant 
properties (see Nesterov and Nemirovsky 1994). Therefore 
from a “warm-start” for a fixed k > 0, i.e. having an approx- 
imation in the Newton area for y(Jt) = argmin {F(y,e,b) | 
y £ R n }, one can improve the current approximation in 
[l-0(m“ 0,5 )] time (m > n) by performing one Newton step 

and updating the barrier parameter k to Jt[l — O(m~ 0 ‘ 5 )] 1 
(see Gonzaga 1992; Nesterov and Nemirovsky 1994). All 
this means that to reach the hot-start, one must perform 
G(v^lnlk 7 ) Newton steps. 

Beginning from the hot-start, we do not need to change 
the barrier parameter k from step to step; instead, we update 
the Lagrange multipliers after at most O^nlnc -1 ) Newton 
steps by minimizing F(y t A, Jb) in y under fixed k > fc 7 . Every 
Lagrange multiplier update leads to an improvement of the 
current approximation by a factor 0 < 7 < 0.5. 

In other words, from the hot-start on, one must perform 
at most 0(ln In ) Newton steps instead of 0(y/rn) Newton 
steps in the interior point methods to improve the current ap- 
proximation by the same factor. Moreover, from some point 
on, the number of Newton steps between two successive La- 
grange multiplier update decreases after every update. This 
is very important for structural optimization because solv- 
ing problem (4)-(5) is only one part of the procedure. The 
second part, which is very time-consuming, is structural anal- 
ysis. This analysis starts with the current approximation for 
y, which we use to recompute the matrix Q and then to solve 
the problem (4)-(5) again. The existence of the hot-start 
makes the optimization part of the procedure substantially 
easier. In the next section, we will show some numerical re- 
sults. 

6 Examples 

The algorithms based on the NMBM to obtain a minimum 
weight structure were used to optimize three truss structures. 
These structures have been used before for studying the char- 
acteristics of the optimality criteria algorithm (see Berke and 
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Khot 1987; Morris 1982; and others). In this investigation, 
the structures were also designed by using the ASTROS pro- 
gram and another program based on an optimality criterion 
(OC) approach. 

In the program using an optimality criterion approach, a 
linear recursive relation based on the reciprocal design vari- 
ables was used to update the design variables, and a linear set 
of simultaneous equations was used to calculate the Lagrange 
multipliers associated with the active constraints. The three 
programs were run on three different computer systems and 
locations and, as such, an exact comparison of the computa- 
tional time was not feasible. The main objective was instead 
to study the convergence characteristics of the methods and 
the various associated programs. 

The analysis package for finite element analysis and eval- 
uation of the sensitivities for the programs based on the op- 
timality criteria and the NMB method were identical. The 
design surfaces near the optimum for all the problems were 
nearly flat and consequently convergence near the optimum 
was very slow. The initial designs for all the structures were 
feasible and were obtained by analysing the structure with 
cross-sectional areas equal to 1.0 in 2 for all the elements, 
then scaling the design to satisfy the constraints. The elas- 
tic modulus was 10 7 lbs/in 2 and the density of the material 
used to calculate the weight of the structure was 0.1 lbs 3 /in 3 . 
The cross-sectional areas of all the elements were treated as 
design variables. 

6.1 Problem 1 - \ 2-element truss 

The structure shown in Fig. 2 was subjected to 4-40 kips at 
nodes 5 and 6, and +20 kips at nodes 11 and 12, respectively, 
in the vertical direction. The displacements at nodes 6 and 
12 were limited to 10 in and 12 in, respectively, and the 
minimum size constraint was 0.1 in 2 . The initial sizes of all 
the elements were 39.746 in 2 and the weight of the structure 
was 57111.67 lbs. The iteration history for the three methods 
is given in Fig. 3. The optimum design weight was 17425.72 
lbs. 



6.2 Problem 2 - 148-element truss 

The truss shown in Fig. 4 was subjected to 40 kips loads ap- 
plied in the x,t/, and z-directions at nodes 7 and 8. Similarly, 
a load of 20 kips was applied at nodes 31 and 32 in the x, y, 
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Fig. 3. Iteration history for 42-element truss 

and z-directions. The displacements at nodes 8 and 32 were 
restricted to 10 in and 20 in, respectively, while the minimum 
size constraint was 0.1 in 2 . The iteration history is given in 
Fig. 5. The initial areas of all the elements were 8.07 in 2 , and 
the weight of the structure was 29474.29 lbs. The optimum 
design weight was 10499.41 lbs. 



Fig. 4. 148-element truss 


6.3 Problem 3 - 721-element truss 

The geometry of the structure is shown Fig. 6. It has 182 
nodes and 721 elements. The length, width and height of 
the box-beam are 600 in, 240 in and 30 in, respectively. The 
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Fig. 5. Iteration history for 148-element truss 


box-beam has six bays in the transverse direction and twelve 
bays in the longitudinal direction. 


cn 



A load of 40000 lbs was applied in the three coordinate 
directions at nodes 25 and 26. At nodes 181 and 182, a load 
of 20000 lbs was applied in the three coordinate directions. 
The minimum size constraint was 1 in'*. 

The initial area of the members was equal to 32.323 in^, 
which was obtained by scaling the structure to satisfy the 
displacement constraints. The initial weight of the structure 
was 114384 lbs. The iteration history Cor this problem is 
shown in Fig. 7. The ASTROS program was terminated after 
twenty iterations because of the excessive computational time 
needed for each additional iteration. The optimum design for 
this structure had a weight 20624.2 lbs. 



Fig. 7. Iteration history for 721-element truss 

The main difference between the three structures was the 
number of design variables. The objective in this study was 
to study how the number of design variables influences the 
convergence characteristics. Comparing the results for ail the 
three problems, it is seen that the initial rate of decrease in 
weight with each iteration was substantially higher for the 
NMBM algorithm than for the other methods. This was 
true for all three problems in spite of the differences in the 
number of design variables. This indicates that the NMB 
method makes maximum use of the information provided by 
the analysis and sensitivity subroutines. It should be re- 
marked that the computational effort for these structural op- 
timization problems, satisfying separability conditions, is the 
least by far for the OC method. The NMBM pays for its 
quick convergence with many internal iterations within an 
outer iteration. ASTROS uses Taylor series approximations 
for internal iterations that are not needed at all for OC. For 
all three methods, each outer iteration is defined here as just 
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one full structural analysis. 

7 Summary 

The Newton modified barrier method, based on rigorous 
mathematical derivation, was used to solve three structural 
optimization test problems. The varying number of size vari- 
ables for the three problms was used to study the effect of 
the number of variables on the number of iterations needed 
to achieve similar levels of convergence. This dependence was 
shown to be weak for the NMBM, similar to this well-known 
characteristic of the basic OC methods. The same problems 
were solved by the ASTROS program with MicroDot as the 
optimizer, and also by an OC-based program. The NMBM 
is shown to have superior convergence properties, reducing 
the objective function at a much faster and more uniform 
rate, thus indicating superior robustness. This fast converg- 
ing robust approach has promising potential for solving struc- 
tural sizing problems with large numbers of variables within 
the emerging multidisciplinary optimization (MDO) problem 
setting. 
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